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Abstract 

It is shown that intrinsic neutrino flavor violation invariably occurs when neutrinos are created 
within the SM augmented by the known massive neutrinos, with mixing and nondegenerate masses. 
The effects are very small but much greater than the naive estimate Am?/E'l or the branching 
ratio of indirect flavor violating processes such as /U — > 67 within the SM. We specifically calculate 
the probability (branching ratio) of pion decay processes with flavor violation, such as vr — > fiv'e, 
showing nonzero results. 
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I. INTRODUCTION 



After the confirmation that neutrinos are massive, nondegenerate and mix themselves, 
further investigations are being intensively carried out, experimentally as well as theoreti- 
cally, to clarify the remaining mysteries about the neutrinos and the new physics they could 
be hiding One question that massive neutrinos immediately poses concerns the status 
of lepton number (L) and family lepton numbers {L^, L^^L^-) that were automatically con- 
served in the standard model (SM) without right-handed singlet neutrinos. We know from 
the successful observation of neutrino flavor oscillations that family lepton numbers are not 
conserved quantities due to the presence of the nondiagonal MNS mixing matrix. Total 
lepton number could be conserved at the classical level if neutrinos were Dirac fermions 
but that scenario does not explain the smallness of neutrino masses. The most natural 
way to explain tiny neutrino masses is the seesaw mechanism but, in this case, neutrinos 
are Majorana fermions in general. Although, approximate lepton number conservation can 
be achieved, guaranteeing small active neutrino masses, by assigning appropriate lepton 
numbers to heavy SM gauge singlets [2I. 

Indeed, it is exactly in the seesaw scenario that many interesting physics could be poten- 
tially observable. If the seesaw scale is relatively low, at the order of TeV, effects such as the 
violation of unitarity of the MNS matrix [3I] may be observable or the direct production of 
heavy seesaw particles 2,|4|, including heavy neutrinos (type I or III) or heavy scalars (type 
II), might be possible. Non-standard interactions could also modify the standard oscillation 
formulas 51]. In such context, it is common to think that all consequences of the SM aug- 
mented by massive neutrinos have been investigated through. (An extensive analysis can be 
found in Ref.O.) Most of the direct consequences of massive neutrinos, with the exception 
of neutrino oscillations, are very difficult to be observed due to the tiny masses and mass 
differences: jAm^gl ~ 8 x 10~^eV^ and |Am23| ~ 2.3 x lO^^eV^ [7]. For example, the pro- 
duction of antineutrinos with negative helicity is possible in principle, because neutrinos are 



massive, but negligible in practice [8]. Despite such difficulties, an enormous experimental 
effort is being dispended to measure the absolute neutrino mass scale jsj. On the other 
hand, indirect effects allowed by massive neutrinos with mixing, such as the lepton flavor 



(LF) violating decay yU — > 67, are even strongly suppressed in t 



because of the tiny neutrino masses that enter the loops 10 



he SM [Br(/i 67) < 10 
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Extensions of the SM, 
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though, may lead in general to relatively large LF violating effects and certain conditions 



should be fulfilled for a natural suppression 



121. 



Contrary to usual expectations, we will show in this article that intrinsic neutrino flavor 
violation, hence lepton flavor violation, is possible in neutrino creation due solely to the 
known neutrino mass differences and nonzero mixing. More specifically, we will show that 
processes such as vr ^ fiUg, are possible with a branching ratio much greater than loop 
induced processes such as /i ^ 67. In fact, this effect should be correctly quantified be- 
fore considering new physics contributions that could mimic the same effects jisl. [l^ . For 
instance, there were attempts to exp lain the LSND anomaly 15|] from new physics inter- 
actions that violate lepton flavor IGj. For interactions that conserve total lepton number, 
however, conflicts with low energy phenomena can not be avoided 17|]. Before the confirma- 
tion that neutrino oscillations were responsible for both deficits of neutrinos coming from 
the sun and the atmosphere, there were attempts to explain the deficit with non-standard 
interactions [Q], even with massless neutrinos 19]. Indeed, it is important to distinguish the 
intrinsic lepton flavor violation effect calculated here from effects coming from interactions, 
extrinsic to the presence of neutrino masses, that violate lepton flavor and, perhaps, lepton 
number. Such interactions could give rise to effective operators with observable consequences 
in other low energy phenomena. An analogous distinction between direct and indirect CP 



20] 



violation is important to classify the CP violating effects involving the neutral i^'-mesons 
that confirmed the CKM mechanism of CP violation in the SM 21]. 

The outline of the article is as follows: in Sec. [IT] we apply the Wigner-Weisskopf approx- 
imation to treat the pion decay, considering the finite decay width. Section HTTl contains the 
main results of neutrino flavor violation in pion decay and uses mainly Eqs. (11 01) and (130|) 
from Secini The ones only interested in the results may skip Sec.[Tll We discuss the results 
and some implications in Sec. lIVI The appendices show some calculations that were omitted 
through the text and some useful material. 



II. WIGNER-WEISSKOPF APPROXIMATION IN PION DECAY 

Consider the pion decay vr" + uj, z = 1, 2 (/i = e,l2 = yu) and j = 1,2,3. The 

detailed description of this decay will be made by applying the Wigner-Weisskopf (WW) 
approximation method [22] . The WW method is essentially an improved method of second 
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order time dependent perturbation theory which can describe the dynamics of decaying and 
decayed states at intermediate times (exponential behavior). 

To calculate the decaying pion state at any time t, within the applicable approximation 
that only liUj states appear as decay states, it suffices to discover the functions ip and x 

where the spin degrees of freedom are omitted and the states {|7r(p)), |/j(q)z/j(k))}, i,j = 1, 2, 
refer to the free states, eigenstates of Hq, normalized as 

(7r(p')|7r(p)) = S\p-p') , 
Uq!)u,{k')\h{q)u,{k)) = 6%ci-q')6%k-k') . (2) 

The expansion ([T]) means we are restricted to the lowest order of perturbation theory. 
The free Hamiltonian is characterized by the free energy of the states with physical masses 

Ho\nip)) = E^ip)\7c{p)) , (3) 
i7o|/.(q)^.(k)) = {EiM + E,^ik))\kici)u,ik)) , (4) 



where Ea{p) = i/p^ + {a = vr, / j, uj), and we will denote M;. = Mj and M^. = rrij. The 
interaction Hamiltonian is given by 

= — y d^xJ^F^x) + counter terms , (5) 

where is the Fermi interaction Lagrangian. 
Considering the total Hamiltonian 

H = Ho + V, (6) 
we can write a Schrodinger-like equation 

i^^^-Ho)Ht))^^ = |rf^P^^^e-^^V(p)) 

+ Aci^k.^^^^^Me-(^.+^^.)*|/,(q)^^.(k)) , (7) 
= V\^{t)) . (8) 
Contraction with the appropriate states yields 

i-V'(p,t) = ^^(P.^) + 5Z yrfVkx.,(q,k;t)(7r(p)|\/(t)|/,(q)z/,(k)) , (9) 



d 



X.,(q,k;t) = j ci3p^(p,t)(/,(q)z/,(k)|\/(t)|7r(p)) , (10) 



where V{t) = e^^'^^Ve and is a counter term. 
From the initial conditions 



V'(p,o) = ^(p), (11) 

Xi,(q,k;0) = 0, (12) 



we can formally solve 



X., (q, k; t) = -I j\t' j d'p ^(p, t') {h{q)u,{k)\V{t')\7c{p)) , (13) 

and obtain 

^^(P,t) = -^^^(P'^) + Jd'p'd\d\ y^rft'(7r(p)|V(t)|/,(q)z/,(k)) 

X (/,(q)z/,-(k)|\/(t')|7r(p'))V^(p',t') . (14) 

This is the key equation for the WW approximation. 

Notice that only momentum conservation holds for the matrix elements, in particular, 

(/,(q)^/,(k)|V^|7r(p)) = iV,;^/'^,,(53(p-q-k) , (15) 

where N^, = (27r)32Ejq)2E,^.(k)2E,(p) and ^i,- = ^i,(p,q,k) = ^(^-(p) ^ Zr(q)i7^.(k)). 
Replacing Eq. (JTSl) into Eq. ( fMl) yields 



d 6M'^ 1 /■* 

-V'(p, t) = -^^^(P, t) - ^^-^ dt'^Pip, t - t')K{p, t') , (16) 



where 

73^ j3 



where Ai?jj = Et, — Ei. —E^_. and the respective p, q, k dependence of Et,, E^, E^. is implicit. 
The expression in Eq. (1171) . however, does not provide a convergent integral since 



behaves as k^ for q = p — k and |k| oo. However, a cutoff function /(p, q, k) multiplying 
is understood to regularize the expression. Such function can arise effectively from the 



pion form factor and vertex corrections in higher orders 231]. Such cutoff function is necessary 
to ensure the convergence of Eq. ( fT7|) and the production rate of vr(p) — > /j(q)z/j(k) to be 
more probable for the energy conserving states and do not grow indefinitely for high |k|. We 
will assume that the cutoff function / satisfies the properties 
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(PI) the functional form of / is broad for Ei. or E^y^ and it varies very slowly for values 
close to the energy conserving values, in particular / = 1 for AEij = 0. 

(P2) the suppression of high momentum |k| or |q| (with q+k fixed) occurs only significantly 
at an scale A which satisfies F ^ A ^ M^/T, where F is the pion decay width. 

Only these properties will be necessary for most of the calculations in this article. The 
inclusion of an explicit cutoff function will be considered in appendix IXlto justify the property 
(P2). 

With the introduction of / we can argue that the dominant contribution of K{p, t) is 
for if: ~ 0, since Eq. (1171) corresponds to a Fourier transform in E^y. and the integrand is a 
very broad function, which leads to a narrow function in time. We can then approximate 
Eq. ([H]) as 



dt' K{p, t') 



^M. (18) 



The Eq. (llSp corresponds to the WW approximation and it is valid for intermediate times, 
i.e., t should be greater than the time width of -ft'(p, t), since for such short time the original 
expression (fT6!) can be significantly different. Within the WW approximation the expression 
inside the bracket in Eq. ffT^ gives 

Using the relation 

^ pi ^ mb{E) , (20) 



E±ie E 

we obtain 



Using the property (PI) of / we can identify Eq. (^T]) as proportional to the pion decay rate 



at rest 



23] 



Re Eq. ^ = M^T , (23) 
while Eq. ( |22l) can be absorbed by the counterterm 

ReEq (dnD = -5M2 . (24) 



We can finally find the functions ip and x- Equation (fT^ gives 

^^(P,^) = -^^(P,^) , (25) 

which can be readily solved to give 

V.(p,t)=^(p)e-^*/2^ (26) 

in accordance to the expected exponential decay law. The factor 7 = Et^{p)/^-k accounts 
for the Lorentz dilatation of time. At the same time, the production wave function can be 
obtained from Eq. (11 01) 



Xii(q,k;t) = x*i(p,q,k;t)^(p)|p^q^j^, (27) 
X,-(p, q, k; t)^[l- e--(Ai..--r/2.).]^-i/2 /^-(p,q, k) _ 



Thus |xij(q, k;t)p is the production probability density. 

From the conservation of probability at any time t, we must check if 

Jd'pmp,t)\' + J2 Jd'qd\\x^MXt)\' = l. (29) 

The calculation is performed in appendix [Bl The important point is that Eq. ( 129|) is satisfied 
if we neglect the terms that does not conserve energy in the squared amplitude I^^P, i.e., 
the second term in 

J2 = l-^f I' + I'^^.jf . (30) 

spins 

where the upperscript EC stands for energy conservation. Notice that the usual energy 
conserving term is positive definite while has no definite sign. The cutoff 

function / is responsible for controlling such contributions. Therefore we retain only the 
energy conserving parts of further on. 

For future use, we also define 

and 

r. = E r,, . (32) 
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The ratio Fj/F corresponds to the branching ratio of the reaction it ^ li + u, independent of 
neutrino flavor, and it practically coincides with the usual branching ratio calculated with 
massless neutrinos, since Ylj = 1 the kinematical contribution of neutrino masses 
are negligible. Obviously, ^ • Fj = F. 

As a last remark, we should emphasize that nowhere in this section the precise form 
of the interaction was used, except in the asymptotic behavior of Therefore, this 

approximation can be used in any two-body decay for which the interaction Hamiltonian is 
known, as long as a proper cutoff function is understood. The explicit amplitude and 
squared amplitude for pion decay are shown in appendix O 



III. NEUTRINO FLAVOR VIOLATION IN PION DECAY 



We begin by defining the two-particle states with definite flavor 24] 



\U(Dufsik)) = Sa^UpJUc^)u,ik)) . (33) 

The charged lepton states remain as mass eigenstates while the neutrino states are mixed 
through Uf3j. We will see, in accordance to usual expectations, that pions decay mainly into 
the states |/a(q)z//3(k)) with (a, (3) = (/x, fi). However, we will also see that there is a non-null 
probability of the pion to decay into the neutrino flavor violating states with {a, (3) = (//, e) 
or {a,j3) = {e,fi). For that purpose, we want to ultimately calculate the probability 

Vi^u.it) = jd\j d'k \Uci)M^)Ht))^j' . (34) 

spins 

Using Xij in Eq. (ITOl) . when t ^ 1/F, we can rewrite Eq. (IMll as 

-Pi^u.it) = ld^ |^(p)p| d'k I J2 Uo.,e-^'''^%F^,f^^^_^ , (35) 



spins J 



where 

TT F f 1 N _ Ar-l/2/'^aj(p,q,k) 

UajFaj (p, q, k) = N^- ^ . (36) 

^ 27 

We see the exponential e~'^'^^* is responsible for the neutrino oscillation phenomenon. In 
fact, if we neglect the neutrino mass mj in every term of Eq. (1351) . except in the exponential, 

we get 

Vi^u,{t) = I d^mp)\'l d'kP,^,^{t)\F^\\ (37) 
8 



where = (-Foj)mj^o- Notice the usual oscillation probability, 



|2 



factors out from the creation probability of \Fa\ , for massless neutrinos. Such factor- 
ization is what allows the definition of the state Eq. (l33l) as a flavor state, since 

Vi^.^it) ^ 6^p^ , (39) 

for 1/r ^ t ^ Lose, where Lose is the typical flavor oscillation length (period). Therefore, 



the antineutrino flavor state Uaj\i^j) is only created jointly with the charged lepton /„ |25l.l26|. 
Notice Eq. ( l39l) correctly coincides with the branching ratio of the decay tt ^ Neutrinos, 
however, are not strictly massless and we may have initial flavor violation because different 
neutrino masses contribute differently to each channel n li + uj [25]. We will focus on 
initial flavor violation and denote the interval of time satisfying l/F ^ t ^ Lose by t = 0. 

We can make the flavor violating contributions explicit by rewriting the term inside the 
square modulus in Eq. (J35l) as 

3 3 
J2 UajU;jFaj = dafsFal + J] UajU^jAF^j , (40) 
j=l 3=2 

where AF^j = F^j — F^i- Thus the square modulus becomes 

3 3 3 2 

l5];[/«,-?7^^.F,,f = 5,/3|F„i|2 + 5«^2Re[F:,^f/,,?7|^.AF„,] + | f/.,f/|, AF„, | . (41) 

j=l j=2 0=2 

We recognize that only the last term of Eq. fl4ip is flavor non-diagonal. The second term, 
which is flavor diagonal, is estimated in appendix [E] and shown to be much smaller than the 
flavor violating contribution. 

Specializing to a 7^ /5, under the approximation of UasU^^ ~ (which is valid if a = e or 
f3 = e), the initial creation probability yields 

'Pi.u,{0) = jd^\^{p)\^jd^k\Uo.,U;,\^\AF^2\'. (42) 

For the two family parametrization, we have |f/a2L'^2p ~ ;| sin^2^, thus indicating that this 
phenomenon is indeed mixing dependent. 

To analyze the most dominant contribution to Eq. (1421) . we recall that a general function 
g{x) can be expanded 

g{x + a) — g{x — a) ^ g'{x)2a , (43) 



for small enough a. Moreover, if g{x) = YYi=i9ii^)y relative difference can be written 



g{x + a) — g{x — a) 

Taking x to be = + 1 Am^ and a = | Am^ 
contributions that compose Fa2'- 

Am^ 



(44) 



|(m2 — ^i) we can estimate the different 



(51) g = E, 



-1/2. 



1^2 



9 
a— 



(52) g = (AE„2 



m 9 



lulEC. 
I^l2 • 



zr/27)-i: 
g' Am? 



9' 
a — 

9 



, Am^ 



7, 



(54) 



2k2 ' 

Am^ /Mf 



m2^0 



2M2 



2M2 V M2 - Ml 



We are assuming the energy conserving values {AEij ^ 0), which is an excellent approx- 
imation considering E^,, |k| are essentially the same either if we compute it using or 



Conventionally we will take the simple average 



^{m1 + m|). In particular, in 



(53), |k|f*^ denotes the momentum of neutrino z/2, assuming energy conservation and p 0: 

(M2-M^+m2)/(2M^). We also note that T > An? /2E^ 



Ik^C 



El^ — ml, where Ejy^ 



is satisfied recalling T = 2.53 x 10~^eV (27| and Am^/2E^ ~ | x 



is either lAm^gl ~ 



10- where Am^ 



0.8 X 10-^ or lAml.l ~ 2.5 x lO'^ [Tj. This condition is necessary to 
have coherent flavor neutrino creation From V <^ E^, |k|^*^, it is also clear that among 
the different contributions (5n), the dominant contribution is given by (52). A through 
analysis of the difference between the amplitudes ^aj, estimated in (54), is shown in ap- 
pendix O Therefore we can neglect all differences due to Am^ in Faj except in the terms 
{AEaj — iT/2j)~^ and obtain 



a2 



AE. 



a2 



N, 



2E,. 



2 ' 



(45) 
(46) 



where the subscript + means we assume m 



m2, as well as in E^ = \/l? + m^. It is 



also implicit that refers to I^^^^P with m] — >• m^ and without the mixing matrix 
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element |f/ojP [see Eq. (|36l) ]. Notice we are already assuming ^ m^, otherwise the term 
inside parenthesis should be kept as E^^(k) — E^-^^(k). Although the |k| — > limit of such 
term in Eq. (1461) is well defined and gives |Am^/ V w? ~ Am = m2 — mi. In practice, for 
realistic |k|, we could assume massless neutrinos for these terms. 

The flavor violating creation probability in Eq. can be calculated in analogy to 
Eq. (lB4p . using Eq. (IHT]) . which gives 

1 r / Am^ \ 2 



^',..,(0)«^s„A.^^— j^^, (47) 

where the two family parametrization, |t/a2f^|2p ~ | sin^26', was employed and p ~ (pion 
at rest) was considered by adjusting ip{j>). The following integral was also necessary, 

i^TT^-^it)- (48) 



One can recognize the term inside parenthesis in Eqs. fjmj) and fHHj) as the additional con- 
tribution that appears in Eq. ( H71) . 

Let us estimate some specific flavor violation probabilities (branching ratios): 

sin22^i2 sin22ei2 T ' sin22^23 

To compute the last value in Eq. (H^ . we considered jAm^gl ^ |Am23| ^ |Am^2l- 



IV. DISCUSSIONS 

The important point of this detailed calculation is that lepton flavor violation should 
necessarily occur when neutrinos are created because it is unlikely that the expression in 
Eq. (H2I) would cancel exactly. It is also important to emphasize that neutrinos should be 
detected as flavor states, as defined (approximately) in Eq. ( l33i) . to observe the flavor vi- 
olation effects. The coherent creation of neutrino flavor states is indeed guaranteed from 
the observations of neutrino oscillations. When neutrinos are not explicitly detected, their 
effects can be computed from an incoherent sum of the contributions of each neutrino mass 
eigenstate 6l, as in the intended direct measurements of absolute neutrino mass. Extensive 

n 

investigations in such context, were first reported in Ref.lH. On the other hand, if mass 
eigenstates were created and detected incoherently, flavor violating effects would be analo- 
gous to flavor changing processes for quarks, at tree level, without the explicit appearance 
of the Am^ dependence. 

11 



The neutrino flavor violation effects reported here are, in general, very small but relatively 
larger than what would be expected from a naive estimate Am^/ (a similar result is indeed 



obtained in Refs. 



28 and 



29l ) because of the presence of the finite decay width F, which is 



very small for pions. We could define, differently from Eq. ( l33l) . that the neutrinos created 
jointly with the charged lepton a is by definition [26!] . However, the difference between 
such definition [Eq. (3.16) of Ref.l26la] and the usual definition in Eq. (1331) carries the factor 
lS.rr? jEl and it is negligible compared to the factor we have calculated in Eq. fHTl) . Thus 



the effect calculated in 
from different sources 



because otherwise there woulc. 
would be no flavor oscillation 



q. WJ\ is dominant, even if we distinguish the neutrinos created 
26|]. In fact, intrinsic neutrino flavor violation effects can not be large 
be no coherent creation of neutrino flavor states and there 



25 



genuine non-standard interactions [l3| . 

Previous calculations of intrinsic neutrino flavor violation 



n 



l26l. 1301 1 . Of course, this analysis is modifled if there are 

r1 



26 



28 



29| did not explicitly 



considered the contribution of the flnite decay width of the parent particle and either ne- 



glected the effect 261] or considered it unphysical 28|, 



29| . The arguments of Ref . l28l is based 



on a formalism that uses a unitarily inequivalent vacuum that guarantees initial neutrino 



flavor conservation 28| but also implies slightly different oscillation formulas 31|- Instead, 
the intrinsic flavor violation effect calculated in Eq. (H7|) should be regarded as a genuine 
physical consequence of massive neutrinos with mixing and it contradicts neither the weak 



Hamiltonian as stated in Ref. 



29 



nor any experimental observations. The qualitative occur- 



rence of intrinsic neutrino flavor violation, that in the context of flavor oscillations could be 



called initial flavor violation, could be anticipated in more phenomenon 
of flavor oscillation probabilities considering scalar 32(] or fermionic 33 



ogical calculations 



34 ] wave packets 



but its magnitude could not be determined without the full consideration of the interaction 
responsible for neutrino creation. 

The expression in Eq. (2) reminds the AS* = 2 contribution from box diagrams in K^- 
mixing (see, e.g.^ Ref.l23l. p. 235). Such contribution is suppressed by the GIM mecha- 
nism [35] because it involves the sum of the contributions of quarks m, c and t in the loop. 
Equation fHTl) . however, is not loop suppressed and, differently for quarks, the mixing angles 
are large. These facts explain the relatively large effect calculated in Eq. (H7|) . despite tiny 
neutrino mass differences. In fact, the effect is much larger than loop suppressed effects such 
as the lepton flavor violating decay /x ^ 67 in the SM. Although, in models beyond the SM 
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such as the MSSM, such effects can be larger than the current experimental limit 111 ]. 

Despite the arbitrariness of the cutoff function /, the expression in Eq. fH2|) is finite 
independently of the presence of that function. This feature shows the robustness of the 
calculation as the cutoff scale A may be chosen from a wide range without affecting the 
results. A expression very similar to Eq. (1471) was estimated in an unrealistic exactly solvable 



QFT model of Lee- type in Ref.|25L also showing that the intrinsic neutrino flavor violation 
effects calculated here bear some universality independently of the particular interaction in 
question. 

The further inclusion of radiative corrections to the formalism developed in Sec.[T] does 
not seem to be straightforward. The corrections have to be included without spoiling the 
conservation of probability of Eq. (1291) . It is also possible that deviations from the exponential 
decay law would emerge from such corrections or from an approximation scheme distinct 
from the WW approximation. Deviations from exponential behavior are indeed expected 
for very short or very long times from the unitary evolution of quantum mechanics 361]. A 
brief connection with perturbative QFT is also shown in appendix [Fl The inclusion of finite 
widths in perturbative QFT is interesting in its own right because it mixes up different 
orders in perturbation theory and special care is necessary in gauge theories to keep track 
of gauge invariance |37l] . Obviously, to fully specify the dominant cutoff scale A, radiative 
corrections should be explicitly considered. The study of the renormalization procedure also 
needs careful analysis. In this respect, it should be emphasized that the necessity of the 
cutoff function / is not related to the nonrenormalizability of the Fermi interaction. The 
same asymptotic behavior ~ k^) would require a cutoff function if instead we adopted 

a Yukawa-type interaction which is renormalizable. In the context of neutrino propagation 



38, 



39| at lowest 



and oscillation, the inclusion of finite widths was also considered in Refs. 
order. 

Another possible application of the formalism developed in Sec. HI] concerns the study of 
the effects of the finite width to the effective size of the decaying particles. The roles played 
by the finite width and the intrinsic momentum uncertainty, encoded here in the wave func- 
tion V'(p); cire not clear but they are crucial to the occurrence of neutrino oscillations, a 
phenomenon that requires quantum coherence. As it is well known, a small uncertainty in 
the spatial localization of the neutrinos are necessary to the observation of neutrino oscilla- 
tions 30|]. With such formalism, the quantum entaglement can be also studied, differently 
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of the static Lee-type model 25j. The extension to three-body decays should be also pur- 
sued since most of the decays with neutrino creation, such as the beta decay or yU — > euei^^, 
have three decay particles. In that respect, it is important to notice that the kinematics of a 
three-body decay is very different from a two-body decay that emits monoenergetic particles 
when the parent particle is at rest. 

To summarize, intrinsic neutrino flavor violation should occur when neutrino flavor states 
are created. The effect is the consequence of the slightly different creation amplitudes, 
functions of different neutrino masses, that have to be summed coherently. The smallness of 
the effect explains why neutrino flavor is an approximately well defined concept in the SM 
and it is directly related to the smallness of the neutrino mass differences. At the same time, 
small mass splittings allow the coherent creation of neutrino flavor states that is required 
for the phenomenon of neutrino flavor oscillations. The observation of the latter enabled 
the recent progress in understanding some of the fundamental properties of neutrinos. 



APPENDIX A: THE CUTOFF FUNCTION / 

We will show here that the contribution to Eq. ( ]B2[) coming from the second term of 
Eq. flHCTj) is negligible for a cutoff function / that obeys the properties (PI) and (P2). We 
will adopt the particular function 

Close to the energy conserving values the contribution of is negligible compared to 

as we can see in Eq. ( ICSp . To analyze the contribution of \6^ij\'^ for |k| ^ F, we 
rewrite Eq. (]C7I) as 

|(5^,,f = {AEi,)^A2 + AE,,Ai , (A2) 

and note that the coefficients A2, Ai are bounded functions of |k|. More specifically, Ei-{p — 
k) — E,^^{k) ^ Ei.{p — k) — |k| is a monotonically decreasing function bounded by -^^^(p) 
and — IpI cos 6*. We also notice that the term inside parenthesis in Eq. ( 1B3I) is bounded as 
well as |1 - e-^(^^«^-^r/27)<|2_ j^us, in analogy to Eq. (IB2]), if we use |/5^ijf instead of 
inside the integral in dk, we recognize we have to compare 



(A3) 
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dX 1,^^ , (A4) 



with 

where A = —AEij and we are neglecting the neutrino masses. Taking only the contribution 
of A2 in Eq. (lA2|l . the ratio between Eqs. (IMIl and (lAil) is 

< |/l2UaxAV 27r7\-i r 

~ l^f P ^ r ; y_[A2 + ^][A2 + A2] ^ ^ 

Ar 

« 1 ' (A6) 



2{M^-Mf 

assuming p and (P2) is valid. The contribution coming from Ai is much smaller 



APPENDIX B: CALCULATION OF EQ. ([29]) 
The second term of Eq. (129 p can be rewritten as 

Jd\£k\xiM,^;t)\' = Jd'pmp)]' Jd\\xiAp.p-K^;t)\' , (Bi) 

where we used the change of variable q — p = q + k and the sum over spins is implicit. 

We can calculate, using (PI), the second integral of Eq. fIBip assuming that the contri- 
bution of the second piece of Eq. (!30l) is negligible and noticing that the squared amplitude, 
summed over spins, after imposing energy conservation, is a function only of the masses: 



2E^{2ny J '\2Ei2E,^dX)^cj_^ p ^2 ^ 

I^^^P /• _ / P dk\ TT r p.. 



recalling that 7 = Et^/M-j^. In Eq. (EH]), the change of variables |k| ^ A = Ey. +Ei. was used 
and the lower end of the integral was extended to —00, considering M^r — Mi — mj ^ F. 

Comparing Eq. (El) with Eqs. (EI]) and (EI]), after using r — ^ijy ^® -^l- (12 9 P is 
satisfied. 
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APPENDIX C: PION DECAY 



The effective Fermi interaction lagrangian is 

££p = -2V2GF{li{x)-i^LUijU,{x))j^{x) + h.c, (CI) 

where L = i(l — 75), {Uij} denotes the MNS matrix while is the hadronic current that 
in the case of pion decay reads 

JfM = VudULlfidi ■ (C2) 
Using we can calculate in Eq. (|T5l) : 



= iCUijUi^qj^^LvjCk) = iCf/jj^jj(p, q, k) , (C3) 



where C = 2Ft,G pVud- We have used the chiral current relation [23 1 



(0|^(x)7,75rf(x)|7r-(p)) = -rV2F^-^M==-^^^ , (C4) 

where ~ 92MeV is the pion decay constant. It is important to keep in mind that 
Eq. (]C3[) should be calculated without assuming energy conservation. In that case, the 
squared amplitude is 

J2 |^.,(p,q,k)| = ^p-qi)ip-k,) - 2(g,-fc,y , (C5) 

spins 

where = (£'^(p),p), gf = (i?;. (q),q) and kj = (i^j,^. (k), k). If we consider energy 
conservation, we get the usual 

J2 \^^AP, q, k)|EC = l^f I = (M,2 - ) + m%M'^ + 2Mf - m]) , (C6) 

spins 

without neglecting the neutrino masses. The remaining part of Eq. (ICSp that does not 
conserve energy can be calculated by using p = q + k + 5p, m. four- vector notation, where 
5p={AE,j,Q): 

\8J^i,\^ = (AE,,)'[P' - {Ei^ - Eu,f] - (M? - m])2AE,,{Ei^ - E,^) . (C7) 
It is important to estimate 



I^Ecp Ml-Mf\M^ Ml) 
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considering A£^jj ~ F, p ^ 0, and the energy conserving values for the rest of the terms. 
Numerically Eq. fICSp is dominated by F/M^ ~ 10^^^ which is negligible and it supports 
why we neglected the contribution of the terms when computing the flavor violation 

probability in Eq. ( H71) . 
We can also calculate 

Aj{p, q, k)^*,(p, q, k) = A{q,-p){p-{k),j,) - 2p\qi-{k)j,>) , (C9) 

spins 

where {k)jj' is given by Eq. ( ]D2I) . To calculate Eq. (1C9I) . we made use of the complete- 
ness relation in Eq. ( IDip . Furthermore, the mixed squared amplitude in Eq. ( ]C9p can be 
decomposed, in analogy to Eq. (I30l) . as 

q, k)^,*,(p, q, k) = \j^f^^,f + P . (CIO) 

spins 

Energy conservation (EC) assumes the neutrino four-momentum is {k)jji. Since the mass 
associated to (/i:)i2, for example, is ^J {k)\2 = 1/^1^2, i-e., the geometrical average, we can 
also show that 

\^l^\^<\^^g\^<\^i^\\ (Cll) 

for mi < Equation (IClip confirms that the contribution due to neutrino mass differences 
in the amplitudes can be indeed neglected in comparison to the contribution containing 
F, i.e., (52), when computing Eq. (1471) . 

APPENDIX D: COMPLETENESS RELATIONS FOR SPINORS WITH DIFFER- 
ENT MASSES 

To compute Eq. (135!) exactly, it is necessary to calculate mixed squared amplitudes such 
as Ylispin ■^ix^iii where the subscripts 1 and 2 denote spinors involving different masses, m\ 
and m^. We are interested, however, in calculating the sum over spins using a common basis 
for the spin directions for the spinors Vy^{\^) and Vy^i^i). (Depending on the parametrization 
adopted Vy^i^^r^Vv^^^s) 7^ 5rs) The only basis where that is possible is their common 
helicity basis. In that basis we have, with helicity 

^t;,,(k,/i)t;,,(k,/i) = (^)2i-7°A^2i, (Dl) 

h 
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where 

{k)2i = y/mim2 ( cosh ^, k sinh ^ , (D2) 
A/C21 = y/mim2 ( cosh A,^, k sinh A^) , (D3) 

with ^ = + ^2) and A^ = ^2 — ^i- The usual hyperbohc parametrization is employed, 
i.e., kj = mj(cosh^j, ksinh^j), with the additional constraint misinh^i = m2sinh^2 = |k|. 
Notice that Eq. ( IDll) reduces to the usual ^ — m when mi = 1712 = m. 
To calculate Eq. (IDll) we made use of the parametrization 



.;,^,(k,/.) = ^^J^t;o(k,/i), (D4) 
where Ej = {kj)^, and the completeness relation 

^t;o(k,%t(k,/.) = l(l-7°). (D5) 

h=± 

APPENDIX E: FLAVOR CONSERVING EFFECTS 

Let us estimate the effect of the second term of Eq. fBT]) which is flavor diagonal. Com- 
paring to Eqs. (H6l) and (H71) . it is possible that it could be relatively large, of the order of 
Am? / EyV . However, we can calculate 

2Re^.;,AF„, « WSa!^^^^-^ . (El) 



After integration in A = —AEa+i the effect is non-null only because of the lower integration 
limit, Aq = mj + Ei^{'p) — Et^{'p) ^ — M^^ + Mq, is finite. One can see the contribution will 
be proportional to 



' r ^ 2E, Ao [A2 + Xi]2 2E,T2tx\1 
which is, in general, much smaller than ( 2e^t ) ^ • 



APPENDIX F: CONNECTION WITH PERTURBATIVE QUANTUM FIELD 
THEORY 

The treatment of unstable states in pertubative Quantum Field Theory (QFT) is of 
considerable interest since the majority of particles studied at high energies, including the 
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ones we call elementary, are unstable and observed as resonances. We will briefly show the 
connection between the formalism developed in Sec. ([11]) with perturbative QFT. 

Although the asymptotic "in" and "out" states can not be defined for an unstable state, 
we can still calculate 

^{p,t) = o(7r(p)|f/(t,0)|7r^)o, (Fl) 
X.,(q,k;t) = o{k{cO,Uj{k)\U{t,0)\n^)o, (F2) 



where 



k^)o = J rf=^pV'(p)k(p))o, (F3) 



and 

ft 



U{t,t') = Texp -i V{t) , (F4) 



with T being the time ordered product. Recall that the S matrix is given by f/(oo, — oo). 
The functions in Eqs. (IF ID and flF2p can be identified with the functions introduced in Sec.HTl 
In particular, they obey the initial conditions of Eqs. (fTTl) and (fT^ . They also obey 

z^V'(p,t) = o(vr(p)mt)[/(t,0)|7r^)o, (F5) 
^^X^,(q,k;t) = o(/^(q),^^i(k)|V(t)t/(t,0)|7r^)o. (F6) 

In particular, if the completeness relation in Fock space could be truncated by 1 = |7r) (7r| + 
we recover Eqs. ([9]) and ( ITOl) . 
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